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Abstract 

Background: Butanol production directly from C0 2 in photosynthetic cyanobacteria is restricted by the high 
toxicity of butanol to the hosts. In previous studies, we have found that a few two-component signal transduction 
systems (TCSTSs) were differentially regulated in Synechocystis sp. PCC 6803 after butanol treatment. 

Results: To explore regulatory mechanisms of butanol tolerance, in this work, by constructing gene knockout 
mutants of the butanol-responsive TCSTS genes and conducting tolerance analysis, we uncovered that an orphan 
slr1037 gene encoding a novel response regulator was involved in butanol tolerance in Synechocystis. Interestingly, 
the Aslr1037 mutant grew similarly to the wild type under several other stress conditions tested, which suggests that 
its regulation on butanol tolerance is specific. Using a quantitative iTRAQ LC-MS/MS proteomics approach coupled 
with real-time reverse transcription PCR, we further determined the possible butanol-tolerance regulon regulated 
by Slrl 037. The results showed that, after slrl 037 deletion, proteins involved in photosynthesis and glycolysis/ 
gluconeogenesis of central metabolic processes, and glutaredoxin, peptide methionine sulfoxide reductase and 
glucosylglycerol-phosphate synthase with stress-responsive functions were down-regulated, suggesting that Slrl 037 
may exhibit regulation to a wide range of cellular functions in combating butanol stress. 

Conclusions: The study provided a proteomic description of the putative butanol-tolerance regulon regulated by 
the slrl 037 gene. As the first signal transduction protein identified directly related to butanol tolerance, response 
regulator Slrl 037 could be a natural candidate for transcriptional engineering to improve butanol tolerance in 
Synechocystis. 
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Background 

Although currently ethanol constitutes 90% of all biofuels, 
biofuels offer a diverse range of promising alternatives. 
Other fuels with better chemical properties, such as buta- 
nol with superior chemical properties in terms of energy 
content, volatility, corrosiveness, and its compatibility with 
the existing fuel storage and distribution infrastructure 
[1], have been proposed as the next-generation biofuels to 
replace gasoline, diesel, and jet fuels. Butanol is typically 
produced by fermentation processes employed anaerobic 
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Gram-positive bacteria, such as Clostridium acetobutyli- 
cum, through an acetone-butanol-ethanol (ABE) fermen- 
tation process [2,3]. However, the relatively slow growth 
rate, complicated life cycle (spore-forming), and significant 
production of by-products associated with the anaerobic 
microorganisms have together made butanol production 
by fermentation processes less competitive with gasoline 
[4,5]. To address these issues, efforts have been recently 
made to produce butanol in more user-friendly non-native 
hosts, such as Escherichia coli and Lactobacillus brevis 
[6-10]. More recently, photosynthetic cyanobacteria have 
attracted a lot of attention for their potential application 
as microbial factories' to produce biofuels and fine chemi- 
cals, due to their natural capability to utilize solar energy 
and C0 2 as the sole energy and carbon sources for growth 
[11]. Lan and Liao (2011) constructed a modified CoA- 
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dependent 1 -butanol production pathway in cyanobac- 
terial Synechococcus elongatus PCC 7942 and achieved 
1 -butanol production from C0 2 and light [12]. By sub- 
stituting a Afunctional aldehyde/alcohol dehydrogenase 
(AdhE2) with separate butyraldehyde dehydrogenase (Bldh) 
and NADPH-dependent alcohol dehydrogenase (YqhD), 
the researchers further increased the C0 2 -derived 1- 
butanol production fourfold [13]. While these exciting 
progresses demonstrated a promising future for renewable 
production of butanol from C0 2 , the current production 
level is still very low, partially due to the toxicity of buta- 
nol towards cyanobacterial cells [14,15]. Therefore, how to 
increase the inherent butanol tolerance of the cyanobac- 
terial production systems may represent an important 
hurdle of improving their productivity [15]. 

Recent genome-level studies have showed that microbes 
tend to employ multiple resistance mechanisms in dealing 
with the stress of a single biofuel product [15-17], and it is 
thus very challenging to achieve tolerance improvement 
by sequential multigene modifications. To address the 
issue, approaches have been proposed to focus on various 
regulatory systems for biofuel tolerance improvement [18]. 
For example, overexpression of the spoOA gene, encoding a 
regulator of stationary-phase events and required for tran- 
scription of solvent formation genes in butanol-producing 
Clostridium acetobutylicum, was shown to enhance the bu- 
tanol tolerance [19]. Random mutation and expression of 
an exogenous global regulator irrE from radiation-resistant 
Deinococcus radiodurans in E. coli resulted in improved 
biofuel tolerances of E. coli cells from several to a hundred 
fold [20]. Random mutagenesis of a global transcription 
factor cyclic AMP receptor protein (CRP) of E. coli resul- 
ted in a twofold growth rate increase when grown under 
1.2% 1-butanol [21]. These studies demonstrated that 
direct manipulation of regulatory systems could provide a 
useful tool in tolerance engineering. However, cur- 
rently very little information is available on regulatory 
systems controlling biofuel tolerance in photosynthetic 
cyanobacteria. 

Two-component signal transduction systems (TCSTSs) 
are the predominant signal transduction system widely 
distributed in bacteria [22]. Interestingly, so far only one 
master response regulator, SpoOA, has been found to be 
involved in butanol tolerance in C. acetobutylicum [19], 
and no TCSTS has been found to be involved in butanol 
tolerance in E. coli or other model species [22,23]. Consid- 
ering that TCSTSs have been shown to be involved 
in sensing and responses to a variety of environmental 
stresses [22], we expect that response to various biofuels 
may also be controlled directly or indirectly by TCSTSs. 
Recently, our laboratory has been focusing on exploring 
tolerance mechanisms of cyanobacterium Synechocystis sp. 
PCC 6803 (hereafter Synechocystis) to various biofuels 
[24-29]. Through a quantitative iTRAQ LC-MS/MS 



proteomics and RNA-seq transcriptomics analyses of 
Synechocystis under butanol stress, we previously identi- 
fied several proteins/genes of TCSTSs differentially 
regulated by butanol exposure [26,29]. In this work, by 
constructing the gene knockout mutant and conducting 
phenotypic analysis, we confirmed that a novel TCSTS 
response regulator Sir 1037 was involved in butanol toler- 
ance in Synechocystis, As the first signal transduction 
protein confirmed related to butanol tolerance, Sir 1037 
could be a potential engineering target of regulatory pro- 
teins for improving butanol tolerance in Synechocystis 
[30]. Towards this end, we further determined the possible 
butanol-tolerance regulon regulated by Slrl037 using a 
quantitative iTRAQ proteomics approach coupled with 
quantitative reverse-transcription PCR (RT-qPCR). 

Results and discussion 

Butanol tolerance analysis of selected TCSTS mutants 

A recent survey showed that Synechocystis contains at 
least 91 genes encoding TCS components, including 27 
histidine kinases, 20 hybrid-type histidine kinases, and 
44 response regulators [31]. Our previous proteomic and 
transcriptomic analysis found four histidine kinases and 
seven response regulators differentially regulated by bu- 
tanol [26,29]. In this study, two TCSTS response regula- 
tor genes, sill 708 and sir 1037, and one transcriptional 
regulator gene, sll0690 } were selected for construction of 
knockout mutants and validation of their involvement in 
butanol tolerance. Currently, no functional information 
is available for sir 1037 and sll0690, while sill 708 (also 
called Rrel7 in the literature) was found to be involved 
in perception of salt stress and transduction of the signal 
in Synechocystis [32]. After chromosomal integration and 
full segregation as confirmed by PCR and sequencing 
(data not shown, the PCR primers for mutant construction 
and validation are listed in Additional file 1: Table SI.), the 
mutants were grown in parallel with the wild-type Syne- 
chocystis in both normal BG11 medium and the BG11 
medium supplemented with 0.25% (v/v) butanol. Com- 
parative analysis showed that there is no visible difference 
in terms of growth patterns between the wild-type strain 
and all three mutants in the normal BG11 medium, which 
was also confirmed using flow cytometric analysis where 
no morphologic difference was detected (data not shown); 
however, gene disruption of sir 103 7 led to an almost two- 
fold increase in butanol sensitivity, suggesting that the gene 
was involved in butanol resistance (Figure 1). The growth 
rates during the exponential growth phase (12 to 36 h) 
were determined as 0.00496, 0.00497, 0.00037, and 0.00250 
OD 630 change per hour for the wild type and the Aslrl037 
mutant grown under normal BG11, and the wild type and 
the Aslrl037 mutant grown under 0.25% butanol, respect- 
ively. No obvious change in terms of butanol tolerance was 
observed for the knockout mutants of sll0690 and sill 708 
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Figure 1 Growth of the wild type and the mutant hslrW37 of 
Synechocystis in BG11 media with or without butanol. BG1 1 

medium was supplemented with 0.25% (v/v) butanol. 



genes (data not shown), suggesting that their differential 
regulation by butanol exposure revealed by previous prote- 
omic and transcriptomic analysis may be part of the 
secondary responses in cells [26,29,32]. 

The sir 1037 gene is also annotated as a CheY subfamily 
RrelO response regulator based on protein domain ana- 
lysis; however, scant functional clues are provided in the 
literature, although previous large-scale protein-protein 
interaction analyses suggested that it could work together 
with a histidine kinase Hikl {sir 1393) [33,34]. Genomic 
context analysis showed that no histidine kinase was 
identified in its neighborhood, suggesting that sir 1037 is 
an orphan response regulator gene. In addition, two genes 
upstream of sir 1037, which were transcribed in the oppos- 
ite direction (s 110985 and s 110984), and three genes down- 
stream of slrl037 (sll0983, sll0982 } and sll0981), which 
were also transcribed in the opposite direction, all encode 
functionally unknown hypothetical genes, providing no 
functional clues for sir 1037. A National Center for 
Biotechnology Information (NCBI) BlastP search using 
the deduced Slrl037 protein sequence showed that the 
genes shared 61%, 58%, 58%, 57%, and 54% identity with 
the response regulators from several other cyanobacteria, 
Cyanothece sp. PCC 7822, Cyanothece sp. ATCC 51142, 
Nostoc sp. PCC 7120, Anabaena variabilis ATCC 29413, 
and Crocosphaera watsonii WH 0003, respectively. How- 
ever, none of the response regulators in these species was 
so far functionally characterized. 

To determine whether the Aslrl037 mutant was sensi- 
tive specifically to butanol, experiments were conducted 
to compare growth patterns of the wild-type strain and 
the Aslrl037 mutant under several other stress conditions: 
1.75% ethanol (v/v), 4% NaCl (w/v), pH 6.5, pH 11.0, and 
8 mM urea. The growth experiment was repeated three 
times, each with three biological replicates. The results 
showed no visible difference between the wild type and 



the Aslrl037 mutant under the tested stress conditions 
(Figure 2). In addition, we calculated the growth rate for 
the strains under various conditions, and the results also 
showed no significant difference between the wild type 
and the Aslrl037 mutant under the tested stress condi- 
tions, suggesting that the sir 1037 gene may not be a 
master response regulator controlling general tolerance 
responses to multiple stresses. It is worth noting that the 
As Ir 1037 mutant behaved similarly to the wild type under 
ethanol stress, implying some degree of the specificity of 
the Slrl037 mediated regulatory network against butanol. 

Overview of proteomic analysis 

For the proteomic analysis, two independent cultivations 
for both wild-type control and the Aslrl037 mutant were 
grown in BG11 medium supplemented with 0.25% butanol, 
and the cells were collected by centrifugation (8,000 x g for 
10 min at 4°C) at the exponential phase (36 h) (Figure 1), 
resulting two biological replicates for the wild-type and the 
Aslrl037 mutant samples. The purified protein samples 
were subjected to the iTRAQ LC-MS/MS proteomic ana- 
lysis as described in the Methods section [26,27]. After data 
filtering to eliminate low-scoring spectra, the qualified 
spectra were matched to 1,521 proteomes, representing 
approximately 43% of the 3,569 predicted proteins in the 
Synechocystis genome [35]. Reproducibility of the quanti- 
tative proteomic analyses were accessed by two types of 
comparisons. First we labeled and mixed two biological 
replicates of either the wild-type or the Aslrl037 mutant 
samples directly for proteomic analysis, and the difference 
was plotted versus the percentage of the proteins identified. 
The results showed approximately 60% of the proteins with 
difference less than a delta error of 0.2 to 0.25, and more 
than 95% of the proteins with difference less than a delta 
error of 0.5 (Additional file 2: Figure S1A,B). Second we 
labeled and mixed each pair of the Aslrl037 mutant and 
the wild-type samples for proteomic analysis, and the 
difference between different biological pairs was plotted in 
Additional file 2: Figure S1C,D. The dispersion of the 
iTRAQ ratios of the quantified proteins was found with 
very similar trends between biological replicates, suggesting 
that the biological noise was reasonably low. 

RT-qPCR validation of proteomic data 

A subset of 13 genes encoding the As/r7037-responsive 
proteins was randomly selected for RT-qPCR validation. 
Among these genes, 6 of their corresponding proteins 
(Slr0179, Slr0476, Slr0879, Ssrl853, S110830, S110108) 
were found to be up-regulated upon the sir 1037 gene 
deletion, and 7 (S110167, Slr0420, Slrl260, S111819, Slr0628, 
S110368, Slrl403) were down-regulated according to the 
proteomic data. RT-qPCR analysis showed a positive 
correlation between RT-qPCR and proteomic results 
for these genes, as 4 out of 6 up-regulated proteins have 
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Figure 2 Growth time courses of the wild type and the hslr1037 mutant of Synechocystis under various stress conditions. A) pH 6.5; 

B) pH 1 1.0; C) 8 mM urea; D) 1.75% ethanol; E) 4.0% NaCI. 



their corresponding genes up-regulated at least 1.5 fold, 
and 6 out of 7 down- regulated proteins have their corre- 
sponding genes down- regulated at least 1.5 fold, respect- 
ively (Table 1). The comparison suggested an overall good 
quality of the proteomic analysis. 

Putative butanol-tolerance regulon mediated by Slr1037 

Using a cutoff of 1.5-fold change plus a statistical signifi- 
cance j?-value less than 0.05, we determined that between 

Table 1 Comparison of ratios calculated from iTRAQ proteomics and RT-PCR analyses 



the Aslrl037 mutant and the wild type grown under 
0.25% butanol, 106 and 164 unique proteins were up- 
regulated and down- regulated, respectively. 

(a) Central metabolic processes 

Our previous transcriptomic and proteomic analyses of 
butanol response in Synechocystis unexpectedly revealed 
that a set of genes/proteins involved in the photosynthesis 
process were up-regulated [26,29]. Similar up-regulation 



Gene ID 


Proteomic ratio* 


RT-PCR ratio** 


Gene description 


slrOI 79 


1.604 ±0.149 


1.197 ±0.140 


Hypothetical protein 


slr0476 


2.335 ±0.791 


1.717 ±0.058 


Hypothetical protein 


slr0879 


1.807 ±0.171 


1 .970 ± 0.098 


Conserved hypothetical protein 


ssr1853 


2.635 ± 1 .097 


1 .857 ±0.1 80 


Hypothetical protein 


SII0830 


1.643 ±0.165 


1 .527 ± 0.094 


Elongation factor g, putative 


sllOI 08 


2.257 ±0.536 


1.131 ±0.152 


Ammonium transporter 


SII0167 


-2.990 ±1.451 


1.210 ±0.1 28 


Hypothetical protein 


slr0420 


-1.838 ±0.504 


-3.448 ±1.101 


Hypothetical protein 


sill 81 9 


-1.811 ±0.402 


-3.980 ± 0.635 


50S ribosomal protein L17 


slr0628 


-1.869 ±0.558 


-1.578 ±0.447 


30S ribosomal protein SI4 


SII0368 


-1.782 ±0.128 


-2.943 ± 0.585 


Pyrimidine operon regulatory protein PyrR 


slr1 260 


-1.975 ±0.539 


-2.303 ±0.570 


Hypothetical protein 


slr1403 


-2.401 ±0.295 


-1.946 ±0.243 


Hypothetical protein 



*Four proteomic measurements were used to calculate the average ratio and standard deviation. 
**Each RT-PCR was repeated three times. 
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of photosynthesis by another biofuel product, ethanol, was 
also found in Synechocystis, and confirmed by increased 
concentration of chlorophyll a concentration in cells [27], 
suggesting that the enhanced photosynthesis related activ- 
ities were related to tolerance to butanol and ethanol in 
Synechocystis. Consistent with this discovery, we found 
that in the butanol-sensitive Aslrl037 mutant, many pro- 
teins involved in photosynthesis activity were down- 
regulated. Six proteins involved in photosystem and 
electron transfer, photosystem II 12 kDa extrinsic protein 
(Sill 194), photosystem II 13 kDa protein homolog (Sir 1739), 
photosystem II reaction center 13 kDa protein (S111398), 
phycobilisome core component (Slrl459), phycobilisome 
small core linker polypeptide (Ssr3383), and cytochrome 
c553 (S111796), were down- regulated (Figure 3). Among 
them, Sill 194 and S111796 were found to be up-regulated 
by butanol, while S111398 was up-regulated by ethanol in 
the wild- type Synechocystis [26,27]. In addition, Slr0399 
was also down- regulated in the Aslrl037 mutant. A 
previous study found that the deletion of slr0399 has 
no significant phenotypic effects other than a decrease in 
thermotolerance under both photoautotrophic and photo- 
mixotrophic conditions, and suggested that Slr0399 is a 
chaperone-like protein that aids in, but is not essential for, 



quinone insertion and protein folding around Q(A) in 
photosystem II, and may also be involved in assembly 
of quinones in other photosynthetic and respiratory 
complexes [36]. 

Two proteins involved in the synthesis of phycobili- 
somes as antenna to enhance light harvesting for photo- 
synthesis in cyanobacteria [37], phycocyanobilin lyase 
alpha subunit (Sir 1878) and phycocyanin alpha phyco- 
cyanobilin lyase CpcE (Slrl098), were down-regulated 
in the butanol-sensitive Aslrl037 mutant. Moreover, 
ferredoxin (S111382) and ferredoximthioredoxin reductase 
(S110554) were also down- regulated. In cyanobacteria, fer- 
redoxins are the main electron shuttles accepting elec- 
trons from photosystem I and delivering them to essential 
oxido-reductive pathways [38], while ferredoximthiore- 
doxin reductase catalyzes the reduction of thioredoxins 
using the [Fe2S2] ferredoxin as a one-electron donor and 
as such plays a central role in light regulation of oxygenic 
photosynthesis [39]. Thioredoxins are small proteins in- 
volved in light-dependent enzyme regulation of many 
oxygen-evolving photosynthetic organisms from cyanobac- 
teria to higher plants [40]. In the Aslrl037 mutant, thiore- 
doxin M (Slr0233) was down-regulated. The protein has 
been found not essential for Synechocystis under normal 
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Figure 3 Schematic representation of central metabolic processes regulated by Slr1037. 
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photoautotrophic growth conditions, but the Aslr0233 
mutant was slightly less tolerant to hydrogen peroxide 
compared to the wild type when grown photoautotrophi- 
cally, and was suggested to have a role in the oxidative 
stress response in Synechocystis [40]. A probable thylakoid 
lumen protein Sill 769 was also down-regulated in the 
As It 103 7 mutant. As a previous study indicated that the 
cytochrome c550 needed to be transported into the thyla- 
koid lumen and then contributed to optimal functional 
stability of photosystem II in Synechocystis [41], the thyla- 
koid lumen protein could be important in maintaining 
proper photosynthesis. 

Two proteins involved in glycolysis/gluconeogenesis, 
phosphoglycerate mutase (Sir 1124) and bisphosphogly- 
cerate-independent phosphoglycerate mutase (Slrl945), 
were down-regulated in the Aslrl037 mutant. The results 
were consistent with those of a previous study, which 
found an increased transcript concentration of phos- 
phoglycerate mutase in cyanobacterium S. elongatus PCC 
7942 during iron limitation [42], and our recent metabolo- 
mic analysis which showed increased 3 -phosphoglycerate 
concentration in Synechocystis under butanol stress [29]. 
These results suggest that an increased withdrawal of 
metabolite flux from the Calvin cycle to increase the rate 
of glycolysis may represent an important mechanism in 
combating butanol toxicity. A protein essential for valine 
metabolism and degradation of leucine and isoleucine, 
3-hydroxyisobutyrate dehydrogenase, was found down- 
regulated in the Aslrl037 mutant (Figure 3). In a recent 
study, 3-hydroxyisobutyrate dehydrogenase was shown to 
confer Pseudomonas putida and E. coli harboring the en- 
coding gene mmsB improved tolerances to the organic 
solvent cyclohexane [43]. 

Compared with the wild-type strain, the As Ir 103 7 
mutant grew poorly in BG11 supplemented with butanol 
(Figure 1). Consistent with a decreased growth pattern, 
proteins involved in core cellular metabolism, such as 
ATP synthesis, protein biosynthesis, protein fate, and 
biosynthesis of purines, pyrimidines, nucleosides and 
nucleotides, amino acids, and sucrose, were among the 
most down-regulated, although most of them are probably 
secondary responses of the sir 1037 knockout rather 
than under direct control by response regulator Sir 1037 
(Figure 3). 

(b) Cellular processes 

An oxidative stress response was observed in E. coli and 
Synechocystis treated with butanol [17,26]. Glutaredoxin 
(Ssr2061) and peptide methionine sulfoxide reductase 
(S111394) were down-regulated in the Aslrl037 mutant 
under butanol stress, suggesting that anti-oxidative 
responses are also regulated by the Slrl037 regulator 
(Table 2). Glutaredoxin has been implicated in mainten- 
ance of a normal cellular thiol/disufide ratio and the 



regeneration of oxidatively damaged proteins [44] and is 
essential for arsenate reduction in Synechocystis [45]. 
Peptide methionine sulfoxide reductase can catalyze 
oxidized methionine (methionine sulfoxides) back to 
methionines, providing cells with a mechanism to repair 
proteins damaged by reactive oxygen species rather than 
having them degraded and then re-synthesizing them de 
novo, S111394 was also up-regulated by ethanol exposure 
in the wild- type Synechocystis [27]. In addition, cyano- 
globin (Slr2097), which has high oxygen affinity, was 
down-regulated in the Aslrl037 mutant under butanol 
stress. A previous study has indicated that cyanoglobin 
sequesters oxygen and presents it to, or is a part of, a 
terminal cytochrome oxidase complex in cyanobac- 
teria [46]. One basic mechanism of salt acclimation in 
Synechocystis involves the accumulation of compatible 
solutes, such as glucosylglycerol [47]. The key enzyme 
in the production of the compatible solute glucosyl- 
glycerol, glucosylglycerol-phosphate synthase (S111566), 
was lower in the Aslrl037 mutant than in the wild type; 
however, the phenotypic analysis showed no change of salt 
tolerance in the Aslrl037 mutant (Figure 2). So far no in- 
formation is available as to whether a similar mechanism 
against salt stress could be used against stress caused by 
butanol or other biofuel products. 

Circadian rhythms of cyanobacteria are controlled by a 
cluster of three genes encoding circadian clock proteins 
{kaiA, kaiB, and kaiC), and the TCS systems involving a 
histidine kinase SasA and two response regulators, RpaA 
and RpaB, have been implicated in their regulation [48]. 
Proteomic analysis showed that KaiB (Slr0757) and KaiC 
(Slr0758) were both down-regulated in the Aslrl037 
mutant under butanol stress (Table 2), implying that 
they are also regulated by the Slrl037 regulator. This 
result is consistent with our previous analysis showing 
that KaiB (Slr0757) was induced by ethanol in the 
wild-type Synechocystis [27]. The potential connection 
between circadian clock and biofuel stress in Synechocystis 
deserves further investigation. 

(c) Cell envelope 

Previous studies revealed that changes of cell membrane 
and envelope represent an important mechanism against 
biofuel toxicity in Zymomonas mobilis, Saccharomyces 
cerevisiae [49], and Synechocystis [26,27]. In the Aslrl037 
mutant, 11 proteins related to cell membrane and enve- 
lope functions were down- regulated (Table 2). Interest- 
ingly, although involved in similar functions, such as 
membrane lipid biosynthesis, protein translocation, and 
formation of the extracellular matrix of the cell mem- 
brane, as those of proteins identified through analyzing 
butanol responses in the wild- type Synechocystis [26], 
down-regulated proteins identified in the As Ir 1037 mutant 
represent a different set of responsive proteins, probably 



Chen et al. Biotechnology for Biofuels 2014, 7:89 
http://www.biotechnologyforbiofuels.eom/content/7/1/89 



Page 7 of 1 3 



Table 2 Proteins down-regulated in the As/r mutant upon butanol exposure 



Protein Mutant_r1 versus Mutant_r2 versus Mutant_r1 versus 
wild-type_r1 wild-type_r1 wild-type_r2 



Mutant_r2 versus 
wild-type_r2 



Description 



Ratio Sequence Ratio Sequence Ratio Sequence Ratio Sequence 
(fold) coverage (%) (fold) coverage (%) (fold) coverage (%) (fold) coverage (%) 



Cell envelope 
SII0274 
SMI 370 
SIM 483 

Slr0015 
Slr0298 

Slr0528 -1.601 21.3 

Slr0774 

Slr0984 

Slrl 271 

SIM510 

Slr1 975 

Slr0528 -1.601 21.3 

Cellular processes 
SII0067 
SII0708 

SIM 394 -1.766 23.9 
SMI 566 

Slr0399 -1.527 40.8 
Slr0757 

Slr0758 -1.65 20.6 

Slr0904 -1.525 10.4 

Slrl 829 -2.171 25.9 



Slr2097 
Ssr2061 



-1J 



55.7 



Other functions 
SII0368 -1.712 60.1 

sin 020 

SII1 249 
SII1 284 
Slr2046 
Ssl2296 

Regulatory functions 
SII1 124 

SMI 672 

SII2012 
Slr0222 



-2.36 1 71 



-1.79 10.9 



-2.16 25.7 



-1.52 25.9 
-1.82 58.6 



-1.65 60.1 



-1.65 11.7 



-1.6 9.4 



-1.52 18.9 
-1.78 8.7 



-1.75 13.8 

-1.79 21.2 

-1.94 11.8 

-1.54 28.7 

-1.99 15.9 



-1.76 23.9 
-1.61 13.8 



-1.95 60.1 

-1.78 18.4 

-1.69 31.9 

-2.01 1.5 



-1.6 



-1.75 



10.3 



-1.541 31.7 



-1.701 21.2 
-1.88 115.3 



-2.237 11.8 



-1.54 15.9 



-1.932 10.9 
-2.464 16.9 



-1.817 60.1 



-1.807 11.7 



-1.761 77.1 



-1.809 2 



HglK-related 

Mannose-1 -phosphate guanylyltransferase 

Periplasmic protein, similar to transforming 
growth factor induced protein 

Lipid-A-disaccharide synthase, putative 

FraH protein homolog 

UDP-N-acetylmuramyl-tripeptide synthetase 

Protein-export membrane protein SecD 

CDP-glucose 4,6-dehydratase (EC 4.2.1.45) 
NAD + dependent 

UDP-N-acetyl-D-mannosaminuronic 
acid transferase 

Fatty acid/phospholipid synthesis protein 
PIsX 

N-acylglucosamine 2-epimerase 
UDP-N-acetylmuramyl-tripeptide synthetase 

Glutathione S-transferase 

Dimethyladenosine transferase 

Peptide methionine sulfoxide reductase 

Glucosylglycerol-phosphate synthase 

Chaperone-like protein for quinone binding 
in photosystem II 

Circadian clock protein KaiB homolog 
Circadian clock protein KaiC homolog 
Competence protein ComM homolog 
Poly(3-hydroxyalkanoate) synthase 
Cyanoglobin 
Glutaredoxin 

Uracil phosphoribosyltransferase 
Probable glycosyltransferase 
Pantothenate synthetase/cytidylate kinase 
Esterase 
MEGF1 

Pterin-4a-carbinolamine dehydratase 



Two-component sensor histidine kinase, 
phytochrome-like protein 

Two-component hybrid sensor and 
regulator 

Group-2 RNA polymerase sigma factor SigD 

Two-component hybrid sensor and 
regulator 
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Table 2 Proteins down-regulated in the As/r mutant upon butanol exposure (Continued) 


SIM 124 


-1.809 2 


Two-component sensor histidine kinase, 






phytochrome-like protein 


Transport and binding proteins 






biiuu64 - 1 .y i \d.o -i.yo 


1 5.8 


Putative polar amino acid transport system 


SII0163 


-1.643 5.7 


WD-repeat protein 


SII0834 -2.36 


11.7 -1.944 11.7 


Low affinity sulfate transporter 


Slr0559 -2.07 


40.7 


Branched chain amino acid ABC transporter 


Slr0864 -1.51 


10.8 


ATP-binding protein of ABC transporter 


Slr0925 


-1.547 35.5 


Single-stranded DNA-binding protein 


Slr1149 -1.734 5 




ATP-binding protein of ABC transporter 


Slr1 302 


-1.786 6.1 


Protein involved in constitutive low-affinity 






C0 2 uptake 


Slrl 540 -1.642 18.6 -1.59 


18.6 


mRNA-binding protein 


Slrl 897 -1.54 


18.9 


Periplasmic sugar-binding protein of 






ABC transporter 


Slrl 974 -1.87 


26.1 


GTP binding protein 


SII0064 -1.91 15.8 -1.96 


15.8 


Putative polar amino acid transport system 



reflecting the complexity of the butanol tolerance network 
utilized by Synechocystis. However, it is worth noting 
that the HglK-related protein (S110274) and FraH protein 
homolog (Slr0298) were down-regulated in the Aslrl037 
mutant. The HglK protein has been found necessary for 
localization of heterocyst-specific glycolipids in Anabaena 
sp. PCC 7120 [50], and the Fra protein influences filament 
integrity, diazotrophy, and localization of septal protein 
SepJ in the heterocyst- forming A nabaena [51]. The mech- 
anism of these proteins against butanol toxicity in the 
unicellular Synechocystis may merit further investigation. 

(d) Transport and binding proteins 

Cross-membrane transporters for small molecules have 
been suggested as one important factor against biofuel 
toxicity in various microbes [40,52]. Our previous prote- 
omic analysis also identified nine putative transporters 
with induced expression level by butanol [26]. In the 
Aslrl037 mutant, eight proteins related to transport and 
binding protein functions were down-regulated (Table 2), 
including a branched chain amino acid ABC trans- 
porter (Slr0559), a periplasmic amino acid-binding protein 
(S110064), a putative polar amino acid transport system 
(S110834), a low-affinity sulfate transporter, a periplasmic 
sugar-binding protein of ABC transporter (Slrl 897), and 
two ATP-binding proteins of ABC transporters (Slr0864, 
Slrl 149). Currently, there is no information available 
regarding substrate specificity of these proteins in Syne- 
chocystis, and their exact physiological function in butanol 
tolerance still needs further investigation. Two systems 
have been identified in Synechocystis to achieve C0 2 - 
concentrating and active C0 2 uptake, one inducible at 
low C0 2 conditions and another constitutive [53]. Our 



analysis showed that Slrl302, a protein involved in cons- 
titutive low-affinity C0 2 uptake, was down-regulated in 
the Aslrl037 mutant (Table 2). Similar regulation of two 
bicarbonate transporters (Slr0040 and Sir 15 12) was also 
found in Synechocystis treated with ethanol [27]. In ad- 
dition, a WD-repeat protein (S110163) was down-regulated 
in the Aslrl037 mutant (Table 2). Synechocystis possesses 
several open reading frames encoding putative WD-repeat 
proteins, one of which, the Hat protein, has been found to 
be involved in the control of a high-affinity transport 
system for inorganic carbon when the cells are grown 
under a limiting concentration of this carbon substrate 
[54]. Although the Aslrl037 mutant was growing poorly 
under butanol stress, no protein involved in the C0 2 fix- 
ation pathway was differentially regulated, consistent with 
the constant expression of the pathway in the wild-type 
Synechocystis treated with butanol [26]; however, the re- 
sults are in contrast with the deceased expression of C0 2 
fixation pathway under stress conditions, such as sulfur 
starvation or low phosphate [55,56], implying that it may 
be C0 2 transport rather than C0 2 fixation that is involved 
in butanol tolerance mechanisms in Synechocystis. 

(e) Regulatory proteins 

Four regulatory proteins were down-regulated in the 
Aslrl037 mutant (Table 2), including a phytochrome-like 
sensor histidine kinase (Sill 124), two hybrid- type sensor 
and regulator proteins (S111672, Slr0222), and an RNA 
polymerase sigma factor SigD (S112012). As sir 1037 is an 
orphan response regulator gene, we expected that its cog- 
nate histidine kinase gene may be standalone and affected 
in the Aslrl 03 7 mutant Examination of three differentially 
regulated TCS genes showed that Sill 124 is an orphan 
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histidine kinase necessary for growth under blue light 
[57], Slr0222 is also an orphan functionally unknown histi- 
dine kinase, and sill 672 is clustered with a downstream 
response regulator Rre2 (slll673). Although group 2 sigma 
factor SigD is nonessential in Synechocystis, it was invol- 
ved in acclimation to salt- or sorbitol-induced osmotic 
stresses [58]. 

(f) Hypothetical proteins 

Sixty-eight hypothetical proteins were down-regulated in 
the Aslrl 037 mutant (Additional file 3: Table S2). Notably, 
S110167 and Sir 1403 were down-regulated in all four bio- 
logical replicates. In addition, several groups of hypothet- 
ical proteins with encoding genes clustered together in the 
chromosome were all down-regulated by the gene knock- 
out, including Sir 1001 -Sir 1002, Slrl259-Slrl260-Slrl261, 
and Slrl275-Slrl276. Using RT-qPCR, we also confirmed 
that slrl403 and slrl260 were down-regulated at tran- 
scriptional level (Table 1). Although currently no func- 
tional information is available for these hypothetical 
proteins, their similar change patterns suggested possible 
involvement in butanol tolerance and may be worthy of 
further investigation. 

Identification of possible Slr1037 binding sites 

To further define the possible regulatory network medi- 
ated by Slrl037, we conducted a promoter DNA-binding 
motif search using 500 bp sequences extracted from the 
upstream regions of all the genes encoding the down- 
regulated proteins in the Aslrl 03 7 mutant using the 
Gibbs Motif Sampler software [59,60]. This analysis showed 
that the top conserved motif identified was a palindrome 
containing 10 total sites with the DNA sequence "GGCG 
ATCGCC" (Figure 4) [61,62]. The motif has several con- 
served positions with information bits greater than 0.5 
[60]. The genes associated with the first motif included 
sll0067 encoding a putative glutathione S- transferase, 
ssr2061 encoding a glutaredoxin 3 (grx3), and sir 1149 
encoding an ABC transporter, which have been sug- 
gested as being related to stress responses in cyano- 
bacteria [44,45,63,64]. As a preliminary prediction, the 
functionality of the motif associated with the Slrl037 
response regulator will need more experimental proof. 



Conclusions 

Increasing studies have found that tolerance to biofuels 
may involve responses of multiple cellular functions 
[14,15]; thus, it may be challenging to achieve biofuel 
tolerance by engineering a single metabolic gene, enzyme, 
or even a pathway. Alternatively, the manipulation of 
regulatory genes could provide a route to complex pheno- 
types, since the regulatory genes tend to control a set of 
genes related to certain functions [18]. To improve biofuel 
tolerance in the photosynthetic cyanobacterial hosts that 
can be used to produce carbon-neutral products directly 
from C0 2 , in this study, we conducted experiments to 
identify regulatory genes related to butanol tolerance in 
Synechocystis. The efforts allowed identification of sir 1037 
encoding an orphan response regulator that was involved 
in butanol tolerance in Synechocystis. Using a quantitative 
iTRAQ LC-MS/MS proteomics approach coupled with 
quantitative RT-qPCR, we further determined the possible 
regulatory network mediated by Slrl037. The results 
showed that along with the sir 1037 deletion, proteins in- 
volved in photosynthesis and glycolysis/gluconeogenesis 
of central metabolic processes, and glutaredoxin, peptide 
methionine sulfoxide reductase, and glucosylglycerol- 
phosphate synthase with stress-responsive functions 
were down- regulated, suggesting that Slrl037 exhibited 
regulation to a wide range of cellular functions in combat- 
ing butanol stress. Finally, it should be noted that, in order 
to cover all possible targets regulated by Slrl037, we used 
a relatively low criterion in determining responsive 
proteins, and the proteins differentially regulated in 
any of the replicates were reported. Thus, more work 
is still needed to further validate whether these molecular 
targets are directly controlled by the Slrl037 response 
regulator. Nevertheless, as the first regulatory protein 
discovered so far that was related to butanol tolerance in 
Synechocystis, the sir 1037 gene can serve as an important 
target for engineering transcriptional machinery in order 
to improve butanol tolerance and productivity in photo- 
synthetic Synechocystis. 

Methods 

Bacterial growth conditions 

Synechocystis sp. PCC 6803 and gene knockout mutants 
were grown in BG11 medium (pH 7.5) under a light 
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Figure 4 Putative regulatory module identified upstream of the genes encoding As/r7037-responsive proteins. The motif is represented 
by a sequence logo generated by the WebLogo software [62]. 
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intensity of approximately 50 [imol photons m" 2 s" 1 in an 
illuminating incubator of 130 rpm at 30°C (HNY-211B 
Illuminating Shaker, Honour, China) [26,27]. Cell density 
was measured at OD 630 on an ELx808 Absorbance Micro- 
plate Reader (BioTek, Winooski, VT, USA). For growth 
and butanol treatment, 10 mL fresh cells at OD 630 of 0.2 
were collected by centrifugation and then inoculated 
into 50 mL of BG11 liquid medium in a 250-mL flask. 
Butanol of analytical purity was purchased from Merck 
(Whitehouse Station, NJ, USA). The growth experiments 
were repeated at least three times to confirm the growth 
patterns. Cells for proteomics analysis were collected by 
centrifugation at 8,000 x g for 10 min at 4°C. 

Construction and analysis of knockout mutants 

A fusion PCR-based method was employed for the con- 
struction of gene knockout fragments [65] . Briefly, for the 
gene target selected, three sets of primers were designed 
to amplify a linear DNA fragment containing the chloram- 
phenicol resistance cassette (amplified from a plasmid 
pACYC184) with two flanking arms of DNA upstream 
and downstream of the target gene. The linear fused PCR 
amplicon was used directly for transformation into Syne- 
chocystis by natural transformation. The chloramphenicol- 
resistant transformants were obtained and passed several 
times on fresh BG11 plates supplemented with 10 [xg/mL 
chloramphenicol to achieve complete chromosome segre- 
gation (confirmed by PCR). Two TCSTS genes, slll708 
and sir 1037, and one transcriptional regulator gene, 
sll0690, that were previously found differentially reg- 
ulated by butanol exposure [26,29], were selected for the 
construction of gene knockout mutants. The successful 
knockout mutants were confirmed by PCR and se- 
quencing analysis. The comparative growth analysis of the 
wild-type Synechocystis and the mutants was performed in 
100-mL flasks each with 10 mL BG11 medium with or 
without 0.25% (v/v) butanol. The cultivation conditions 
were the same as described above. Growth analysis was 
performed in biological triplicates. 

Proteomics analysis 

i) Protein preparation and digestion 

For each sample, 10 mg of cells were frozen by liquid 
nitrogen immediately after centrifugation and washed 
with phosphate buffer (pH 7.2). The cells were broken 
with sonication cracking at low temperature. The cell pel- 
lets were then resuspended in a lysis buffer (8 M urea, 4% 
CHAPS, 40 mM Tris-HCl) with 1 mM PMSF and 2 mM 
ethylenediaminetetraacetic acid (EDTA) (final concentra- 
tion). After 5 min of vigorous vortexing, dithiothreitol 
(DTT) was also added to a final concentration of 10 mM. 
After mixing, the samples were centrifuged for 20 min at 
20,000 x g, and the supernatant was mixed well with ice- 
cold acetone (1:4, v/v) with 30 mM DTT. After repeating 



this step twice, the supernatants were combined and pre- 
cipitated at -20°C overnight, and stored at -80°C prior to 
sample cleanup if not for immediate use. For digestion, 
protein pellets from the previous step were resuspended 
in digestion buffer (100 mM triethylammonium bicarbon- 
ate, TEAB, 0.05% wlv sodium dodecyl sulfate, SDS) to a 
final concentration of 1 mg/mL (total protein measured 
by bicinchoninic acid assay (Sigma, St. Louis, MO, USA)). 
Equal aliquots (500 \ig) from each lysate were then 
digested with trypsin overnight at 37°C (Sigma; 1:40 w/w 
added at 0 and 2 h) and lyophilized. 

ii) iTRAQ labeling 

The iTRAQ labeling of peptide samples derived from the 
control and butanol-treated samples was performed using 
the iTRAQ Reagent 8-plex kit (Applied Biosystems, Foster 
City, CA) according to the manufacturers protocol. Four 
samples (two biological replicates for the wild-type control 
and two biological replicates of the Aslrl037 mutant, 
respectively) were iTRAQ labeled. The peptides were 
labeled with respective isobaric tags, incubated for 2 h, 
and vacuum centrifuged to dryness. The labeled wild-type 
control and Aslrl037 mutant replicate samples were 1:1 
pooled, and generated four combinations of samples, 
which were reconstituted in buffer A (10 mM KH 2 P0 4 , 
25% acetonitrile, pH 2.85). The iTRAQ labeled peptides 
were fractionated using a PolySULFOETHYL ATM SCX 
column (200 x 4.6 mm, 5 \xm particle size, 200 A° pore 
size) by an HPLC system (Shimadzu, Japan) at flow rate 
1.0 mL min" 1 . The 50-min HPLC gradient consisted of 
100% buffer A (10 mM KH 2 P0 4 , 25% acetonitrile, pH 2.85) 
for 5 min; 0 to 20% buffer B (10 mM KH 2 P0 4 , 25% ACN, 
500 mM KCL, pH 3.0) for 15 min; 20 to 40% buffer B for 
10 min; 40 to 100% buffer B for 5 min followed by 100% 
buffer A for 10 min. The chromatograms were recorded at 
218 nm. The collected fractions were desalted with 
Sep-Pak @ Vac CI 8 cartridges (Waters, Milford, MA), 
concentrated to dryness using a vacuum centrifuge, and 
reconstituted in 0.1% formic acid for LC-MS/MS analysis. 

Hi) LC-MS/MS proteomic analysis 

The mass spectroscopy analysis was performed using an 
AB SCIEX TripleTOF™ 5600 mass spectrometer (AB 
SCIEX, Framingham, MA, USA), coupled with an online 
micro flow HPLC system (Shimadzu, Japan) as described 
before [26,27]. The peptides were separated using a Nano- 
bore C18 column with a PicoFrit nanospray tip (75 \im 
ID x 15 cm, 5 \im particles) (New Objectives, Woburn, 
MA). The separation was performed at a constant flow 
rate of 20 \iL min" 1 , with a splitter to get an effective flow 
rate of 0.2 \iL min" 1 . The mass spectrometer data was 
acquired in the positive ion mode, with a selected mass 
range of 300 to 2,000 m/z. Peptides with +2 to +4 charge 
states were selected for MS/MS. The three most abundantly 
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charged peptides above a count threshold were selected 
for MS/MS and dynamically excluded for 30 s with ±30 
mDa mass tolerance. Smart information-dependent acqui- 
sition (IDA) was activated with automatic collision energy 
and automatic MS/MS accumulation. The fragment inten- 
sity multiplier was set to 20, and the maximum accumula- 
tion time was 2 s. The peak areas of the iTRAQ reporter 
ions reflect the relative abundance of the proteins in the 
samples. For peptide identification, the Triple TOF 5600 
mass spectrometer used in this study has high mass 
accuracy (less than 2 ppm). Other identification parame- 
ters used included: fragment mass tolerance: ±0.1 Da; 
mass values: monoisotopic; variable modifications: Gin- > 
pyro-Glu (N-term Q), oxidation (M), iTRAQ8plex (Y); 
peptide mass tolerance: 0.05 Da; max missed cleavages: 1; 
fixed modifications: carbamidomethyl (C), iTRAQ8plex 
(N-term), iTRAQ8plex (K); other parameters: default. 

iv) Proteomic data analysis 

The MS data were processed using Proteome Discoverer 
software (Version 1.2.0.208) (Thermo Scientific) to gen- 
erate the peak list. The default parameters of Proteome 
Discoverer software (Version 1.2.0.208) were used. The 
data acquisition was performed with Analyst QS 2.0 soft- 
ware (Applied Biosystems/MDS SCIEX). Protein identifi- 
cation and quantification were performed using Mascot 
2.3.02 (Matrix Science, London, UK) [26,27]. For iTRAQ 
quantification, the peptide for quantification was automat- 
ically selected by the algorithm to calculate the reporter 
peak area, error factor (EF), and j^-value (default parame- 
ters in the Mascot software package). The resulting data 
set was auto bias-corrected to eliminate any variations 
imparted due to the unequal mixing during combining of 
the different labeled samples. The genome sequence and 
annotation information of Synechocystis sp. PCC 6803 
were downloaded from NCBI, the Comprehensive Micro- 
bial Resource (CMR) of TIGR, and CyanoBase [35]. Pro- 
teins with a 1.5 fold or more change between the Aslrl037 
mutant and the wild-type control samples and a ^-value 
of statistical evaluation less than 0.05 were determined 
as differentially expressed proteins. The quantitation was 
performed at the peptide level by following the procedures 
described previously [26,27]. The Students t-test was per- 
formed using the Mascot 2.3.02 software. 

Quantitative real-time reverse tanscription-PCR analysis 

The identical cell samples used for proteomics analysis 
were used for RT-qPCR analysis. cDNAs were synthesized 
using RevertAidTM Reverse Transcriptase (Fermentas, 
Glen Burnie, MD, USA). The qPCR reaction was carried 
out in 20 uL reactions containing 10 uL of SYBR® Green 
PCR Master Mix (Applied Biosystems, Foster City, CA), 
and 2 uL of each PCR primer at 2 mM, employing the 
StepOnePlus™ Real-Time PCR System (Applied Biosystems, 



Foster City, CA), under the following conditions: 50°C for 
2 min and 95°C for 10 min, followed by 40 cycles of 95°C 
for 15 s and 60°C for 1 min. Quantification of gene expres- 
sion was determined according to the standard process of 
RT-qPCR, which used serial dilutions of known concen- 
trations of chromosome DNA as a template to make a 
standard curve. The rnpB gene (6803 sOl) encoding 
RNase P subunit B was used as an internal control [66]. 
Three technical replicates were performed for each gene. 
Data analysis was carried out using the StepOnePlus 
analytical software (Applied Biosystems, Foster City, CA). 
The data was presented as ratios of the amount of 
normalized transcript in the mutant to that from the 
wild-type control. The gene ID and their related primer 
sequences used for real-time RT-qPCR analysis are listed 
in Additional file 1: Table SI. 

Promoter analysis and motif identification 

The Gibbs Motif Sampler software from the Biometrics 
Laboratory of Wadsworth Center, Albany, NY, USA was 
used to identify matrix models describing DNA sequence 
motifs present upstream of genes responsive to butanol 
treatment [59,60]. Regions representing approximately 
500 base pairs of the DNA sequences upstream of the 
translational start site of responsive genes to butanol stress 
were extracted from the NCBI genome database using the 
Regulatory Sequence Analysis Tools (RSAT) [61]. Both 
strands of each sequence were searched and possible motif 
locations were identified using the motif matrix score ob- 
tained from the Gibbs Motif Sampler software. The multi- 
level consensus sequence for each motif was then used to 
generate a sequence logo that is a graphical representation 
of nucleic acid multiple sequence alignment [62], 

Additional files 



Additional file 1: Table SI. Primers used in this study. 

Additional file 2: Figure SI. Repeatability between biological 
replicates. Wild-type (A) and I±shl037 mutant (B), respectively. Distribution 
of iTRAQ log ratios of the detected proteins between As/r/037 biological 
replicate 1 and the wild type (C) and between kslr1037 biological 
replicate 2 and the wild type (D). 

Additional file 3: Table S2. Hypothetical proteins differentially regulated 
in the Aslr1037 mutant. 
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